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The grand-canonical potential of many-body physics can be considered as a functional of the 
interaction-free correlation function. As such it obeys a nonlinear functional differential equation 
which can be turned into a recursion relation. This is solved graphically order by order in the two- 
particle interaction to find all connected vacuum diagrams with their proper weights. As a special 
case, the procedure is applied to generate the Hugenholtz diagrams for a weakly interacting Bose 
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PACS numbers: 24.10.Cn 



I. INTRODUCTION 



A systematic and physical approach to construct all Feynman diagrams of a quantum field theory together with 
their weights was proposed some time ago by Kleinert [1]. It is based on a functional differential equation for the 
vacuum energy which involves functional derivatives with respect to free propagators and interactions. Solving the 
_ ■ functional differential equation by a recursive graphical procedure leads to all vacuum diagrams. In a subsequent 
\Q ' step, the diagrams of n-point functions are found graphically by amputating lines and vertices in the vacuum 
diagrams. Recently, this approach was used to systematically generate all connected and one-particle irreducible 
Feynman diagrams of the euclidean multicomponent </> 4 -theory both in the disordered symmetric phase [2] and in 
the ordered, spontaneously broken-symmetry phase [3, 4]. The approach was also applied to QED [5] and scalar 
QED [6]. Furthermore, a modification of this program led to a closure of the infinite hierarchy of Schwinger-Dyson 
equations for the n-point functions in a certain functional sense. Both in 4 -theory [7, 8] and in QED [9] this allowed 
> us to construct most directly all diagrams for n-point functions which are relevant for the renormalization of the theory. 



Here we shall use this approach to generate iteratively the vacuum diagrams of many-body theory. We start with 
briefly reviewing some relevant graphical operations in Section II. They are applied in Section III to determine the 
grand-canonical potential as a power series of the two-particle interaction. In Section IV we derive the functional 
differential equation for the grand-canonical potential and show how to solve it graphically order by order. Finally, 
Section V discusses as a special case the resulting recursive graphical construction of the Hugenholtz diagrams for a 
weakly interacting Bose gas. 



II. MANY-BODY PHYSICS 

The thermal fluctuations of a quantum many-particle system are controlled by the action functional 

1-1./,*./, 1 1 f T/( int )„ 



AWM= / G^V^ + r / V^Vi«V>4, (1) 

J12 L J1234 

where ip*^ denote complex fields for bosons and Grassmann fields for fermions. In this short-hand notation, the 
number indices of the fields ip*,ip, the bilocal kernel G~ x , and the two-particle interaction ]/( int ) are abbreviations for 
the spatial coordinates and the imaginary time, i.e., 

i = {x-i,n}, J^ = J d 3 x 1 J dn, i/)* = ^(xrjTi) , ^ = V(xi,ri) , 

G£ =G- 1 (x 1) t 1 ;x 2 ,t 2 ), = int >( Xl) n; x 2 , r 2 ; x 3 , r 3 ; x 4 , r 4 ) . (2) 



Both for bosons and for fermions, the action functional (1) is specified by the kernel 



G (xi,ti;x 2 ,t 2 ) = 5(xi -x 2 )5(n -r 2 ) 



^-2^ A2 + y(X2) ^ 



(3) 
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where M denotes the mass of the particles and V represents the one-particle potential. Furthermore, the two-particle 
interaction is a functional tensor with the symmetry 

v l23i — v 3412 • W 

In this paper we shall leave the kernel G _1 in the action functional (1) completely general, and insert the physical 
value (3) only at the end. By doing so, we regard the action functional (1) as well as all local and global statistical 
quantities derived from (1) as functional of the bilocal kernel G _1 . In particular, we are interested in studying the 
partition function 

Z = J ViP* J Vip e ~ A ^''^ , (5) 

where the functional integral is performed with respect to those fields ip*,ip which are periodic or antiperiodic in 
the imaginary time t, depending on whether we treat bosons or fermions. By expanding the functional integral 
(5) in powers of the two-particle interaction y( mt ), the expansion coefficients of the generating functional consist of 
interaction-free expectation values. These are evaluated with the help of Wick's rule as a sum of Feynman integrals, 
which are pictured as diagrams constructed from lines and vertices. The interaction-free correlation function G, which 
is the functional inverse of the kernel G _1 in the action functional (1) 



Gi 2 G^ 1 = S 13 , (6) 

J2 

is represented by a straight line with an arrow 

G12 = 1 — 2 , (7) 
and the two-particle interaction y( int ) is pictured as two vertices which are connected by a dashed line 

-« - ><• (8) 

The graphical elements (7) and (8) are combined by an integral which graphically corresponds to attaching a line to 
a vertex as, for instance, 



. >-< S . (9) 

1 4 4 — — 5 1 5 



In this paper we generate the subset of connected Feynman diagrams contributing to the generating functional (5) 
together with their weights. To this end we introduce the functional derivative with respect to the interaction-free 
correlation function G according to 

^=S 13 6 i2 . (10) 
OG34 

Such functional derivatives are represented graphically by removing a line of a Feynman diagram in all possible ways 
[1-9]. For practical purposes it is convenient to use also functional derivatives with respect to the bilocal kernel G _1 
whose basic rule reads 

S M = 6is6 a . (11) 

As has been elaborated in detail in Rcfs. [1-9], such a functional derivative is represented by a graphical operation 
which cuts a line of a Feynman diagram in all possible ways. Indeed, from the identity (6) and the functional product 
rule, we find the effect of this derivative on the interaction-free correlation function 



SG12 



G13G42 , (12) 



which has the graphical representation 

<5 



<5G 3 4 



2 = 1 — - — 3 4 — - — 2 . (13) 
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The functional derivatives with respect to the correlation function G and the kernel G 1 are related via the functional 
chain rule 

■J=T = - f G 3iG 24 Jr- . (14) 
oG 12 J 3i oG 3i 

Thus cutting a line is equivalent to amputating a line and adding two lines to the vertices to which the original line 
was connected. 



III. PERTURBATION THEORY 



Field theoretic perturbation expressions are usually derived by introducing external currents j*,j into the action 
functional (1) which are linearly coupled to the fields i(>*,ip [10-15]. Thus the partition function (5) becomes in the 
presence of the generating functional Z[j* 7 j] which allows us to find all interaction-free n-point functions from 
functional derivatives with respect to these external currents To calculate these, it is possible to substitute 

two functional derivatives with respect to the currents j*,j by one functional derivative with respect to the kernel 
G _1 . This reduces the number of functional derivatives in each order of perturbation theory by one half and has the 
additional advantage that the introduction of the currents j*,j becomes superfluous. 



A. Current Approach 

Recall briefly the standard perturbative treatment, in which the action functional (1) is artificially extended by a 
source term 

A[r,W\j]=AP,4>]- Jiiiih + tth) , (15) 

where denote complex current fields for bosons and Grassmann current fields for fermions. The functional integral 
for the generating functional 

Z\j*,j] = Jvtp* J Thl> e-^'-^'jl (16) 
is first explicitly calculated for a vanishing two-particle interaction V^ lnt \ yielding 

ZW\j*,j] = cxp (^TrlnG- 1 + G 12 f l]2 ^ , (17) 
where the trace of the logarithm of the kernel is defined by the series [16, p. 16] 

00 (— b)™+ 1 r 

TrlnG- 1 = £ { —^- / (G^ S 12 ) ■ ■ ■ (G^ 1 - <5„i) . (18) 

n=l U Jl - n 



If the two-particle interaction l/( mt ) does not vanish, one expands the generating functional Z[j*,j] in powers of 
l/( lnt ), and reexpresses the resulting powers of the fields within the functional integral (16) as functional derivatives 
with respect to the currents The original partition function (5) can thus be obtained from the interaction-free 

generating functional (17) by the formula 



z = cxp I" \ I ^4 , jfrTjA z( ° )[f ' j] 

I 2 J 1234 oj 1 dj 2 dj 3 dj 4 J 

Expanding the exponential in a power series, we arrive at the perturbation expansion 

f y(int) ^ l 1 f^Y f ^(int) (int) 

_ l + ~T / 1234 Ti TvT« JT* + 9 \ ~o~ I / 1234 5678 

[ 2 J 123i 0jidj 2 dj3<3j 4 2 \ 2 J .712345678 



(19) 

j«=0,j=0 



(20) 

j*=0,j=0 

in which the pth order contribution for the partition function requires the evaluation of 4p functional derivatives with 
respect to the currents 
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B. Kernel Approach 

The derivation of the perturbation expansion simplifies, if we abandon introducing the currents j*,j and use instead 
functional derivatives with respect to the kernel G -1 in the action functional (1). In this way we obtain for the partition 
function (5) instead of the previous expression (19) the result 



exp 



{ o / ^1234 c n -l c n -l } e ° ' ( 21 ) 
I 1 J1234 oG 12 dG M J 



where the zeroth order of the grand-canonical potential fl^ = In Z^ has the diagrammatic representation 

q(0) = ipTrlnG" 1 = T Q • (22) 
Expanding again the exponential in a power series, we obtain 

X2 1 / 1\2 



7 - J i + zl / T/ (int) ^ + I fzT\ /" T/ (int V ( 
_ 1 / 1234 Ar^Ar^ 1 2 I IT J / 1234 5< 

z J 1234 0(_r 12 0(_r 3 4 Z \ Z / J 12345678 



int) I fi(°) 



5678 r/^-ix/^-i 
r 12 uv - z 34 " \ ~ / "'12345678 6Lr 12 " ( - T 34 0t *56 6Lr 7& 



e" 1 ' . (23) 



Note the two advantages of this expansion over the conventional one (20) in terms of currents coupled linearly to 
the fields. First, it contains only half as many functional derivatives. Second, in case of fermions it does not contain 
derivatives with respect to Grassmann fields. Taking into account (11), (12), and (18), we obtain 

6 e n(W = TG 21 e- (0) , (24) 



S 2 



SG 12 SG 3 l 



e nm = (G 21 G43±G 23 G 4 i) e o<0) , (25) 



such that the partition function Z becomes 

Z = (l - \ I (G 21 G 43 ± G 23 G 41 ) + . . . 1 

I z J1234 J 

According to the Feynman rules (7) and (8), this is represented by the diagrams 

*={i4o-o 4e +-}«p{to } 



(26) 



(27) 



Note that each closed loop causes a factor ±1. 

IV. GRAPHICAL RECURSION RELATION FOR CONNECTED VACUUM DIAGRAMS 

In this section, we derive a functional differential equation for the grand-canonical potential Vl whose solution leads 
to a graphical recursion relation for all connected vacuum diagrams. 

A. Functional Differential Equation 

The functional differential equation for the grand-canonical potential Vt is derived from the following functional 
integral identity 

J Vr j^-^ (V> 2 * e-WM) = (28) 

with the action (1). This identity is the functional generalization of the trivial integral identity dx f'(x) = for 
functions f(x) which vanish at infinity. Nontrivial consequences of Eq. (28) are obtained by performing the functional 
derivative in the integrand which yields 

J Vr J Vi, (s 12 =F J Gr 3 V 2 >3 T f ^MV-sV^s) e~ A ^^ = . (29) 
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Substituting the field product ^2^3 by a functional derivative with respect to the kernel G 23 according to 

V> 2 > 3 e-WM = --|=t e~ A ^'^ , (30) 

dLr 23 

this equation can be expressed in terms of the partition function (5): 

ZS 12 ± [cti™=±[ C'™- (3D 

J3 0(^23 ^345 OCr 23 0Cr 45 

Going over from the partition function Z to the grand-canonical potential 

n = lnZ, (32) 
the linear functional differential equation (31) turns into a nonlinear one: 

J 3 <$G 23 J345 V <5G 23 £G 45 <5G 23 SG 45 J 

If the interaction potential y( int ) vanishes, this is immediately solved by (22) due to (24). For a non-vanishing 
interaction potential V^ lnt \ the right-hand side in (33) produces corrections to (22) which we shall denote with r2^ lnt ). 
Thus the grand-canonical potential decomposes according to 

n = O (0) + ft (int) . (34) 

Inserting this into (33) and taking into account (24) as well as (25), we obtain the following functional differential 
equation for the interaction part of the grand-canonical potential f^ lnt ); 

r r (int) r (int) f (int) sn(^ 

± / G 13 77^1 ^ / ^1345 ^34^52 ± / V 1345 G32G54 - / ^1345 G 32 ! 
J3 0G 23 7 3 4 5 7 3 4 5 7345 0U 45 

/■ 1/( int )rr 5QP*>,[ ..(int) / ^^ (!nt) ^ ^ (lnt) 

/ l/ 1345 (j 54 ! ± / l/ 1345 I , , ± j x I ■ (OOJ 

J345 0G 23 J345 V 0<-r 23 <)G 45 (5G 23 0G 45 / 

Setting 1 = 2 and performing the integration over 1, we obtain because of the symmetry (4) of the interaction potential 

y(int). 



/■ r -l^ (int) _ f „(int) r r _ /■ „(int) r p ±2 / y (int) r 

/ G 12 - _ / "l234 G21G43 =F / V 1234 G41G23 ± ^ / * / 1234 ^21 -r-~r 

'12 "G 12 J1234 J1234 J1234 "G 34 



r(iat) s 2 n(^ r (int) sn^ sn^ 



^1234 / ^1234 sn -l xn -l ■ (36) 

4 oG 12 0G 34 J1234 d G 12 oCr 34 

Equation (36) contains functional derivatives with respect to the kernel G -1 which are equivalent to cutting lines 
in the vacuum diagrams. For practical purposes, however, it will be more convenient to work with derivatives with 
respect to the correlation functions G which remove lines. With the help of the functional chain rule, the first and 
second derivatives with respect to the kernel G -1 are rewritten as (14) and 



■ c^-l — / G7iG28G5 3 G4 6 — — h / (G51G23G46 + G4lG26G5 3 ) — — , (37) 

0G 12 <3G 34 J5678 OG7 8 OG 56 J 56 0tx 56 



respectively. Taking into account again the symmetry (4) of the interaction potential V^ nt ^ , the functional differential 
equation (36) for fi( int ) takes the final form 



L 



int 



G12 — — — — I ^1234 G21G43 =F / ^1234 ^2-iG'n 



12 0(J12 J 1234 J 1234 



v ant) r r r sn(^) r (int) snw 

^1234 ^51^23^46— ^ T ^ / V1234 "2ll*53W6 — — ; 

0Lx56 ./1934Kfi OG5 



123456 0"~r56 J 123456 0*^56 

<5 2 fi (int) /■ „<mt^ „ „ „ <5^ (int) sn^ 



12345678 



(int) S 2 Q {lnt > f ( int ) £Qt lnt J (5f]( mt ) 

^1234 G7lG28G 53 G4 6 — — / V 1234 G71 G28 G53 G46 -r— ; — — . (38) 

OLr780(jr 56 J12345678 0<^78 OCt 56 
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Supplementing the Feynman rules (7) and (8) with a graphical representation of the interaction part of the grand- 
canonical potential 

tt (int) = (^) , (39) 
the functional differential equation (38) can be depicted graphically as follows: 



c;H-o--o±e+»cH"0-< £ 



1 s 2 



3 <5l^^20~3^^4 (5l^^2 ■> 
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Recursion Relation 



,X.^r- (40) 



We now convert the functional differential equation (40) into a recursion relation by expanding the interaction part 
of the grand-canonical potential into a power series of the interaction potential y( int ); 

oo 

= E(Z> ( 41 ) 

p=l 

The vacuum diagrams which contain p times the interaction potential y( int ) have the property that they satisfy the 
eigenvalue problem 



As the operator on the left-hand side removes a line from the Feynman diagram which is later on restored, it counts 
the number 2p of lines in a vacuum diagram. Inserting the expansion (41) into the functional differential equation 
(40), we obtain for p > 1 the recursion relation 



1 / 1 



. v - y ± 2 r>-< ^ 



2(p+l) I 



and the initial diagrams 



+ : 2 ^+V^ >--<T T^l (43) 

'^T 3 (51^2(53^4 ^ oT^2 ^-4 (53^4 1 

4 9=1 



5O--0 4e ■ («> 



The right-hand side of (43) contains four graphical operations. The first three are linear and involve one or two 
line amputations of the previous perturbative order. The fourth operation is nonlinear and mixes two different line 
amputations of lower orders. To demonstrate the working of this recursion formula, we construct the connected 
vacuum diagrams in second and third order. We start with the amputation of one or two lines in the first order (44): 



2^ W ± - J,< ' 0W2OW4 3^1 ' ^ 



5i^2 2' v_y 2 ^]/ ' (5i^2<53^4 3 ^: 3 ^>; 
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Inserting (45) into (43), where we have to take care of connecting only legs with the same label, we find the second-order 
correction of the grand-canonical potential f2: 

© =3 a::o + o-o ± \ o-o-o ± 10 ± 30- (46) 

The calculation of the third-order correction f^ 3 ) leads to the following 20 diagrams: 

© = 5 g:d + 1 •:>•::• + \ o-o + o::o ± ,' g^d 
+5 o-© + 0-0 + 0-0 + 0-0 ± a-o::o 

4 0-0-0 ± 0-0-0 ± 0-0-0 + 5 qS-o 

40-0-0-0 4© o©±© 4© 40- <«» 

From the vacuum diagrams (44), (46), and (47), we observe a simple mnemonic rule for the weights of the connected 
vacuum diagrams which is similar to a corresponding one for QED [5, 9]. At least up to four loops, each weight is equal 
to the reciprocal number of lines, whose amputation leads to the same two-point diagram. The sign is given by (±)', 
where / denotes the number of loops. Let us also note that the sum of all weights of the connected vacuum diagrams in 
the loop order under consideration vanishes. These simple weights are a consequence of the bose and fermi statistics 
as well as the interaction (8). The weights of the vacuum diagrams in other theories, like <^ 4 -theory [2-4, 17, 18], 
follow more complicated rules. 



V. HUGENHOLTZ DIAGRAMS 

Now we consider the special case of a dilute bose gas where the two-particle interaction VO nt ) can be approximated 
by the local one 

C = y«Wi4. (48) 

Here A denotes the thermal wave length and a the scattering length. This specification modifies the graphical 
representation (8) of the two-particle interaction. The dashed line is removed and the two vertices are united to a 
single one: 

2 >-< 3 => • (49) 

1 S X 4 j / \ 4 

The resulting vacuum diagrams are called Hugenholtz diagrams [10, 19]. In this way the direct and the exchange 
vacuum diagram of the first order (44) of the grand-canonical potential reduce to 



= OO ( 5 °) 



which is also shown as diagram # 1.1 in Tabic I. Correspondingly, the second and the third order (46) and (47) lead 
to the vacuum diagrams depicted in Table I. We observe that the weight of a Hugenholtz diagram for the weakly 
interacting bose gas is given by the formula 

^=2WJV (51) 
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(compare with the corresponding formula for the (/> 4 -theory [2, 17, 18]). Here D,T,F denote the number of the 
following double, triple and fourfold connections between two vertices: 



D 



■■XX' T - • F: ^ • (52) 



Furthermore, the number TV stands for the number of vertex permutations leaving the vacuum diagram unchanged, 
where the vertices remain attached to the lines emerging from them in the same way as before. The Hugenholtz 
diagrams also follow from a graphical recursion relation. From (43) and (49) we get for p > 1: 



2(p+l) I ^^-2 5l^2 ^ V4 61^283 



p- 1 5{ 1 ) 1 ^- 4 6 (p-q) 



<5l^^2 2-^ V -3 S3 
9=1 



This is iterated starting from (50). Indeed, amputating one or two lines in the first order (50), i.e. 



6(1) 1 _ <5' 



= 2 



we obtain from (53) for p = 1 the second-order diagrams in Table I. A subsequent amputation of one line 

5C2 



(Si ^—2 




2 2 



•0- 1 (55) 



and two lines 



5l^^2 8z^^a 



s2(2) ?0 O? iO Oi 

= 4 2 . ^ 1 +4 2 V . / . ■ 1 +4 2 . W 3 +4 4 V . / . . 3 

4 44 2 

leads with p = 2 in Eq. (53) to the third-order contribution in Table I which also shows the results for the subsequent 
order. We remark that the grand-canonical potential for the dilute bose gas has been calculated so far up to three 
loops in the seminal paper [20]. Recently, this result has been used to determine the critical temperature for the 
homogeneous dilute bose gas up to second order [21]. 
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TABLE I: Hugenholtz diagrams and their weights for the grand-canonical potential of a weakly interacting bose gas. Each 
diagram is characterized by the vector (D,T,F; N) whose components lead to the weight according to Eq. (51). 



